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ABSTRACT 

We determine the orientations of the light distribution of individual elliptical 
galaxies by combining the profiles of photometric data from the literature with 
triaxial models. The orientation is given by a Bayesian probability distribution. 
The likelihood of obtaining the data from a model is a function of the parameters 
describing the intrinsic shape and the orientation. Integrating the likelihood over 
the shape parameters, we obtain the estimates of the orientation. We find that 
the position angle difference between the two suitably chosen points from the 
profiles of the photometric data plays a key role in constraining the orientation 
of the galaxy. We apply the methodology to a sample of ten galaxies. The 
alignment of the intrinsic principle axes of the NGC 3379, 4486 and NGC 5638 
are studied. 

Subject headings: galaxies : photometry - galaxies : orientation. 
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1. Introduction 

Intrinsic shapes and orientations of the individual elliptical galaxies have been 
investigated by Binney (1985), Tenjes et al. (1993), Statler (1994a,b), Bak and Statler 
(2000), Statler (2001) and Statler et al. (2004). These authors have used the kinematical 
data and the photometric data, and have used the triaxial models with the density 
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distribution p(m 2 ), where m 2 = x 2 + \ + % with constant axial ratios p and q. It was 
shown analytically that the projected density of such a distribution p(m 2 ) is stratified 
on similar and co-aligned ellipses (Stark 1977, Binney 1985). Statler (1994a,b, hereafter 
S94a, S94b) uses (apart from the kinematial data) a constant value of ellipticity, which 
is an average over a suitably chosen range of the radial distance, for the shape and the 
orientation estimates. The shape estimates are robust, and are described by a pair of shape 
parameters, namely the short to long axial ratio of the light distribution and the triaxiality 
of the mass distribution. The orientations are described by the Euler angles which define 
the orientation of the galaxies relative to the line of sight. 

The methodology is based on Bayesian statistics, and is described in S94a, S94b, and 
more recently in Statler et al. (2004). The likelihood of obtaining the data from a model in 
a chosen orientation is a function of the orientation parameters and the shape parameters. 
The likelihood is multiplied by a model for the parent distribution of these parameters 
(the prior), and integrated over the parameters that one is not trying to constrain. The 
normalized result is the posterior probability distribution for the parameters of interest. 
Probability of the shape is obtained by integrating over the uninteresting orientation 
parameters, and likewise, probability of the orientation is obtained by integrating over 
the shape parameters. The shape and orientation thus evaluated do not constitute a 
simultaneous fit to the data ; they are, in fact, mutually exclusive (S94a). 

Chakraborty et al. (2008, hereafter C08) have demonstrated that the variation in the 
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intrinsic shapes of the light distribution of elliptical galaxies can be determined by using the 
photometric data alone. They use the models with varying axial ratios (Chakraborty 2004) 
and the triaxial models discussed by deZeeuw and Carollo (1996). These models exhibit 
ellipticity variation and position angle twist. They use ellipticities 6i n , e ou t and the position 
angle difference Q out — 0m at two suitably chosen points i? m and R out from the profiles of 
the photometric data of the galaxies. The intrinsic shape is described by the parameters 
(go, Qoo, To, Too), where go, ?oo(To, T^,) are short to long axial ratios (triaxialities) at small 
and at large radii, respectively. They find that the best constrained shapes parameters are 
the axial ratios go and g^, and the absolute value of the triaxial difference T^ defined as 
\Td\ = {Too — Tq\. The methodology adopted in C08 is essentially the same as described in 
S94a with necessary modifications to suit the requirements. 

We now take up the complementary problem. We determine a statistical estimate of 
the orientation of the light distribution of the galaxies using photometric data. We consider 
the usual polar coordinates (6 ,0) of the line of sight, with respect to the principle axes 
of the galaxy, to define the orientation. We use ellipticities e m , e out and the position angle 
difference = Q ou t — 6m, as in C08, and calculate the likelihood of obtaining the data 
from a model in a chosen orientation. 

We consider a flat parent distribution (go, goo, To, Too, 4> ) = constant (flat prior) and 
multiply it with the likelihood to obtain the posterior density. Integrating the posterior 
density over the shape parameters, we obtain the marginal posterior density V(9',(J)) 
describing the statistical estimate of the orientation. It is necessary that the marginal 
posterior density should be likelihood dominated, and not prior dominated. The prior, 
which is considered as flat in the present investigation, may be replaced by an improved or 
modified prior when the results of a large sample of galaxies are available. Modified prior is 
estimated and used in Bak and Statler (2000). 
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Although it is the posterior density which constitutes the result of the Bayesian 
estimates, some statistical summary is also useful for its description. Area of the parameter 
space enclosing 68% of the total probability, which can be interpreted as la error bar is 
one such useful quantity. A working definition of a good likelihood dominated estimate is 
that the 2a region (area enclosing 95% of the total probability) should be less than 50% of 
the total parameter space (S94a). Another useful quantity is the location (9 ,(j) p ) of the 
maximum probability. 

In the present investigation, we consider the triaxial galaxies under study to have very 
simple photometric structure. Each one of them is considered as a triaxial distribution 
of one component luminous matter. Thus, we ignore complications, such as dust, shells 
or disks etc. In this context, we mention that dust lanes or disks are sometimes used to 
constrain the orientations and shapes (see, eg. Tenjes et al. 1993, S94a). We also assume 
that the axes of each triaxial galaxy are intrinsically aligned, and therefore, the orientation 
as determined by using the observed data at a small Ri n and at a large R out gives the 
orientation of the entire galaxy. We shall examine this assumption in section 5. 

Determination of the intrinsic shape and the orientation using photometry is important 
because the number of galaxies with good photometry is many more than those with good 
kinematics. We find that photometry can constrain the orientation of the light distribution 
of the triaxial galaxies. We find that the position angle difference plays a key role in 
constraining the orientation. The most probable line of sight of a galaxy with very small 
|Gd| lies close to one of the principle planes, although the estimates are not very much 
likelihood dominated. The orientation estimates are likelihood dominated and are well 
constrained for galaxies with large |0d|, typically around 4°.0 or more. The marginal 
posterior density is symmetrical in 4 pairs of octants in all cases of small and large |0<j| 
(These points are described in detail in sections 3 and 4). 
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Section 2 presents the triaxial mass models. Section 3 presents the orientations of NGC 
7619 and two synthetic galaxies, which establish our methodology. Section 4 discusses the 
orientations of a sample of galaxies. In section 5, we examine the alignment of the intrinsic 
principle axes of NGC 3379, 4486 and 5638. Results and a discussion are presented in 
section 6. 



2. Model 

We use models which are described in C08. The models are triaxial generalizations of 
the spherical 7-models (0 < 7 < 3) of Dehnen (1993). These models have cusp at the centre 
and the cumulative mass converges. The projected surface density corresponding to 7 = 1.5 
most closely resembles to the de Vaucouleurs law, and we concentrate to 7 = 1.5 models 
only. However, it is realized that the shape and orientation estimates are not very sensitive 
to the radial profile of the projected density (S94b). 

A triaxial generalization of Dehnen's spherical density distribution p(r) is made triaxial 
by considering the distribution p(r) and replacing r by M, where M 2 = x 2 + + ?p- Here, 
(x, y, z) are the usual Cartesian coordinates, r is the spherical radial coordinate and 

{ } [5b 2 + M 2 ' {) 
Q -2 (M) = fl>V + M 2 q~J 

W 1 j 13b 2 + M 2 • [) 

The constant density surfaces are coaxial and coaligned ellipsoids with varying axial 
ratios (P,Q), which are (po,Qo) at small radii and are (pocQoo) at large radii. In above, b 
is the scale length and (3 > is a parameter. The parameter (3 alters the values of (P, Q) 
in the intermediate region. These models are presented in Chakraborty (2004) and are 
modified in C08. We shall refer to these as M 2 models. 
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Another form of triaxial generalization of Dehnen's model is presented in deZeeuw and 
Carollo (1996), and is modified in C08. Rewriting p(r) as f(r), the triaxial model is density 
distribution of the form 

P = f(r) - W) + 9l (r)]Y 2 ° + [h(r) + ^i(r)]F 2 2 , (3) 

where p is the density in the usual spherical polar coordinates (r,6,<j)), g(r) and h(r) are 
two suitably chosen radial functions, Y" 2 ° = 3/2 cos 2 9 — 1/2 and Y" 2 2 = 3 sin 2 9 cos 20 are the 
usual spherical harmonics and 

9i = a . ( Lr _l , 7 [30rr 2 ], (4) 
4ir{r + r 2 J ' 

hi = a ^ [30rr 4 ]. (5) 
47r(r + r 4 )' 

In above, g(r), /i(r) and the parameters r\..r^ are the same as in deZeeuw and Carollo. 
gi{r) and h\(r) are introduced in C08, which for choices of parameter a > 0, alter the 
profiles of e and in the intermediate region. L is the luminosity of the model. We shall 
refer to these as fgh models. 

Both the M 2 and the fgh models form a family of four parameters models, and are 
fixed by a choice of axial ratios (p , qo,Poo, <?oo)- The triaxialility T and are related to 
the axial ratios at small and at large radii by 



1 - Po nr 1 - P A 

^0 = 

for the fgh models, and by 



To = , ^ = —Zf. (6) 



oc 



2 



rp _ 1 PO rp _ 1 Poo 



for the M 2 models. To fix up the scale length b of the triaxial models, we consider the 
value of effective radius R e = 1.286 of the spherical models. The effective radius of the 
triaxial models depends on the axial ratios, as well as on the viewing angles. However, 
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such changes are small for 7 models (deZeeuw and Carollo, 1996) and are neglected. For 
a choice of (go, To, q^, T^), we take various values of /3 and a to calculate the probability. 
The unweighted sum of the probability using these ensemble of models gives a model 
independent estimate of the shape and orientation. 

We study the variation of the projected properties of the models with orientation. We 
consider (1) a M 2 model with = 1 and (g , T , q^, T^) = (0.9,0.1,0.5,0.2) and (2) a fgh 
model with a = and (go, T , g^, T^) = (0.7,0.3,0.9,0.1) and calculate the 'photometric' 
data at Ri n = 0.25R e and R ou t = 1.40i? e (R e = 1.28). To project these models we take a 
fixed value of 9' and a varying </>'. The variations in 6j n , e out and as a function of <p' are 
presented in Table 1. The variation in the position angles difference is noteworthy. The 
variation is large, and it changes sign across <$' = 90°, 180°, 270° and 360°. We find that 
this has a significant role in constraining the orientation. On the other hand, changes in 
ellipticities with 0' is relatively small. 

To gain some further insight in the kind of results that we expect, we plot the 
constant G^ contours as functions of (9',(p'), for models with chosen parameters. Figure 1, 
presents such constant G^ contours for fgh models. Figure 1A is plotted with parameters 
(g , To, goo, Too) = (0.85,0.2,0.70,0.3), while figure IB is drawn for model parameters 
(go, T , goo, Too) = (0.85,0.5,0.70,0.2). Eye inspection of these plots (and the similar plots 
for other choices of the model parameters) indicates the possible range of values of (9' , <$') 
where the observed G^ of a galaxy is close to the calculated G^ from a model. In this range, 
we shall obtain higher values of the likelihood. In particular, we find that the contours of 
small Qd are close to one of the principle planes. 
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Table 1: Variation of projected properties of the models with orientation. 

M 2 Model with (q Q , T , q^, T^) fgh Model with (q Q , T , q^, T^) 

= (0.9, 0.1, 0.5, 0.2) k6' = 30° = (0.7, 0.3, 0.9, 0.1) & 9' = 60° 



rh' 










tout 




i n 

1U 


n nzL 


0, 


.04 


1 3 ^A 
io.o^ 


n m 

U.Ul 





,01 


-1.11 


9(1 
zu 


n nzL 

U.U^ 


0, 


06 


1 A AO 


n m 

U.Ul 





,01 


1 ^3 
-l.OO 




D DzL 


0, 


.08 


1 1 D7 


n m 

U.Ul 





,02 


1 33 
-l.OO 


-tv 


n 

u.uo 


0, 


10 


7 Q9 


n D9 





.02 


-u. yu 


OU 


n 

u.uo 


0. 


12 


7 90 


n D9 





.03 


D fizL 

-U.U^l 


ou 


u.uu 


n 
u 


i 3 
10 




n n3 
u.uo 


n 
u 


03 

,uo 


n zLi 

-U.^l 


70 


0.06 


0, 


.14 


3.60 


0.03 





,03 


-0.24 


SO 


0.06 


0, 


.14 


1.02 


0.03 





,04 


-0.11 


90 


0.06 


0, 


.14 


0.01 


0.03 





.04 


0.00 


100 


0.06 





,14 


-0.37 


0.03 





,04 


0.11 


110 


0.06 





.14 


-2.96 


0.03 





.03 


0.24 


120 


0.06 





.13 


-4.60 


0.03 





.03 


0.41 


130 


0.05 





.12 


-6.49 


0.02 





.03 


0.64 


140 


0.05 





,10 


-7.85 


0.02 





,02 


0.96 


150 


0.04 





,08 


-10.80 


0.01 





,02 


1.33 


160 


0.04 





,06 


-13.10 


0.01 





,01 


1.53 


170 


0.03 





,04 


-12.86 


0.01 





.01 


1.11 



- 10 - 



3. Method 

We apply the methodology outlined in section 1 to determine the orientations of NGC 
7619 and two synthetic galaxies. 

3.1. NGC 7619 

The observed data of NGC 7619 is e in = 0.28 e out = 0.19, 6 d = 4°.0 at R m = 10.55 
arcsec and R ou t = 45.71 arcsec. Effective radius R e of NGC 7619 is 32.0 arcsec. This data 
is obtained from R-band photometry of Franx et al. (1989). While applying this data to 
calculate the probability of orientation, we consider the uncertainty of 0.02 in ellipticities, 
and the uncertainty of 1°.0 in position angles both at R{ n and -Rout- 
Figure 2 presents the plot of the posterior probability distribution V of the orientation 
as a function of polar coordinates (6 , cp') of the line of sight. We use a rectangular plot 
with 6' going from 0° to 180° in the vertical direction from bottom to top and (j) going 
from 0° to 360° in the horizontal direction from left to right. We label the eight octants of 
the space of the parameters (#',</>), for some convenience in discussions. These labels are 
shown in figure 2. The origin of the coordinate system is 'placed' at the centre and the axes 
are 'drawn' along the principle axes of the triaxial galaxy. Here, we have used fgh models 
with a = and have used 6j n , e ou t and B^. The probability is shown in dark-grey shade : 
darker is the shade, higher is the probability. The white contour encloses 68% of the total 
probability which can be interpreted as la error bar. 

It is seen that the higher probability regions are disjoint blobs, each confined to one 
of the octants of the parameter space. Further, the plot is quite symmetric. Part of the 
plot in the pair of octants (1, 2) is repeated in the pairs (3, 4), (5, 6) and (7, 8). We thus 
conclude that the orientation is well constrained and it suffices to consider the probability 





90 o 



90 



Fig. 1. — Contours of constant for fgh model with a = 0. Figure 1A is drawn with 
model parameters (go, T , q^, T^) = (0.85, 0.2, 0.70, 0.3), while IB is drawn with parameters 
(0.85,0.5,0.70,0.2). 




Fig. 2. — Plot of the posterior probability distribution V cis du function of (0 , ) for NGC 
7619. Labels 1 to 8 denote the octants (see, text for detail). 
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max=0.235W9E-02 




Fig. 3. — Plot of the posterior probability distribution V function of (9 , ) for NGC 
7619 in octants (1, 2). Plus marks the location of the maximum probability 




Fig. 4. — Plot of the posterior probability distribution V function of (0 , <p ) for NGC 
7619 using e in and e out only. 
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distribution V in one of these pairs of octants only. 

In figure 3, we present the plot of probability V in octants (1, 2). Here, we use M 2 
models with (3 = 0.2, 1.0, and 5.0 and fgh models with a = 0.0, 2.5 and 5.0 and take 
the unweighted sum of the probability over all these models. We find that the la region 
occupies 27.9% of the total parameter space. The probability is likelihood dominated and 
is well constrained. The location of the maximum probability is at 9 = 40°.5, <j> = 106°.5. 
In addition to above, the la region extends to 3 more pairs of octants, as described above 
and exhibited through figure 2. 

In obtaining the plots of figures 2 and 3, we have used e» n , e out and 9^. We find that 
position angle difference 0^ has a very important role. To illustrate this, we present the 
plot of V of NGC 7619 in figure 4, wherein we use 6j n and e out only. We find that ellipticities 
alone do not constrain the orientation. In particular, almost all values of 4> are in the higher 
probability region. In the orientation estimates presented hereafter, we use the complete 
set of variables, namely 6j n , e out and O^. The position angle difference 0^ plays a major 
role in obtaining a likelihood dominated orientation estimates whereas the ellipticities ej n 
and e out and the choices of Ri n and R out produce a fine tunning of the posterior probability 
distribution. 



3.2. Synthetic galaxies 

In section 3.1, we presented the orientation estimates of NGC 7619. The 
reliability of the results can be tested by comparing them with results of the other 
workers. Alternatively, one can test of the results by recovering the orientation of an 
observed model from the ensemble. We consider the M 2 (/3 = 1,7 = 1.5) model with 
parameters q = 0.9, T = 0.1, = 0.5, = 0.2 and project it at viewing angles 
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9' = 30°. and <jj = 45°. 0. We refer to this as galaxy A. The "observed" parameter 
are e in = 0.0408, e out = 0.0899, O in = 11°. 6, and Q out = 17° .4 at R in = 0.25R e and 
Rout — lA0R e (R e = 1.28). Here 0^ is 5°. 8. The uncertainties in ellipticities and in position 
angles are taken as the same as in the previous case of NGC 7619. Using ei n ,e ou t and 0^ 
the plot of V over the entire range of (9 ,$)■ It is qualitatively similar to figure 2, and 
therefore, similar to the case of NGC 7619, we conclude that it is sufficient to consider the 
probability in one of the pairs of octants. 

Figure 5 presents the plot of V to galaxy A in octants (1, 2). The contour encloses 
the region of 68% of the total of probability. The hedge marks the true orientation and the 
plus marks the location of the highest probability. We find a good agreement between the 
results of probability calculations and the true orientation. This can be regarded as a proof 
of the reliability of the methodology. 

As another example of synthetic galaxy, we consider M 2 model with the same values 
of (go, Tq, goo, Too) as those for galaxy A, but orientated at 9 = 40°. and 4> = 150°. 0. The 
"observed" parameters are e in = 0.0573, e out = 0.1117, Q in = —6°. 4 and Q out = —10°. 9. R in , 
R out and R e are the same as in galaxy A. The position angle difference 0^ is —4°. 5. We 
refer to this as galaxy B. The plot of the probability V{9 ' , ) of galaxy B in octants 1 and 
2, is presented in figure 6. Here again we find that the true orientation (marked by hedge) 
and the location of maximum probability (marked by plus) are close to each other. 

The results of the galaxies A and B indicate that the orientation estimates of the 
galaxies using photometry alone are reliable. We find that galaxies A and B, as well as 
NGC 7619 exhibit large position angle difference (|©d| ~ 4°.0 or more) and the orientations 
of these galaxies are well constrained. 
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Fig. 6. — Same as Fig. 3, but for the synthetic galaxy B. The plus marks the location of 
highest V, while the hedge marks the true orientation. 
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4. Orientations of a sample of elliptical galaxies 

We now apply the methodology, which is now matured enough through its application 
to NGC 7619 and synthetic galaxies A and B, to a sample of elliptical galaxies. We find 
it convenient to divide the sample in categories based on the values of the position angle 
difference. This is based on large, moderate and small values of |0<z|, without any sharp 
boundaries between them. NGC 7619 and the synthetic galaxies presented in section 3 are 
placed in large |0^| category. Table 2 presents the list of elliptical galaxies of our sample, 
alongwith the observed data used in the models. Here L, M and S denote large, moderate 
and small | ©a] categories. 

We begin by examining the orientations of the galaxies with large |0<f|. Galaxies other 
than NGC 7619, having large \G d \ are NGC 1407, 2986, 4374, 4486 and 5638. We present 
the plots of V of NGC 1407, 2986, 4374, 4486 and 5638 in octants (1,2) in figures 7 - 11. 
We find that the orientations of galaxies are well constrained. We present a summary of the 
orientations of these galaxies in Table 3. 

Next, we consider galaxies with small |0<f|. Galaxies NGC 3379, 4261 and 4551 are 
placed in this category. Plot of the posterior density of NGC 3379 in the entire range of 
(9 , (f)') is presented in figure 12. Although it is not constrained (almost all value of (j) are in 
la regions), but the symmetry of the plots in pairs of octants is clearly seen. Therefore, it is 
sufficient to consider the probability in octants 1 and 2 only. Plots of the probability of the 
orientation of galaxies NGC 3379, 4261 and 4551 in octants (1, 2) are presented in figures 
13, 14 and 15 respectively. We find that the orientation is not constrained, specifically, 
almost all values of 4> are allowed in the high probability region. 

Galaxy NGC 1052 with 0^ = 2°. 6 is a border line case between the large and the small 
1 0d | categories. Probability density of the orientation of NGC 1052 is presented in figure 
16. It is seen that the orientation is little poorly constrained : 68% highest probability 
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Table 2: Observational data used in models. 



Galaxy 


R e 










Qd 


category 


Source* 


NGC 1052 


36.5 


10.2 


35.4 


0.269 


0.330 


2.6 


M 


1 


NGC 1407 


72.0 


19.6 


75.1 


0.050 


0.060 


-5.0 


L 


2 


NGC 2986 


41.0 


7.3 


36.4 


0.160 


0.130 


5.0 


L 


2 


NGC 3379 


37.5 


15.7 


49.3 


0.078 


0.133 


0.0 


S 


1 


NGC 4261 


42.5 


10.6 


58.8 


0.251 


0.154 


0.3 


S 


1 


NGC 4374 


57.0 


7.3 


63.6 


0.193 


0.076 


-5.3 


L 


1 


NGC 4486 


110.0 


15.0 


69.0 


0.033 


0.100 


-13.3 


L 


1 


NGC 4551 


21.0 


5.3 


32.4 


0.227 


0.284 


1.0 


S 


1 


NGC 5638 


29.5 


8.0 


36.9 


0.056 


0.118 


8.4 


L 


1 


NGC 7619 


32.0 


10.5 


45.7 


0.280 


0.190 


4.0 


L 


2 



a (l) Pelctier et al. (1990), (2) Franx et al. (1989). 



Table 3: Summary of the orientation estimates. 



Galaxy 






Percentage of the area of 
68% highest probability 


NGC 1052 


79.5 


139.5 


34.1% 


NGC 1407 


13.5 


127.5 


33.0% 


NGC 2986 


31.5 


118.5 


27.4% 


NGC 4374 


67.5 


175.5 


26.5% 


NGC 4486 


73.5 


10.5 


11.0% 


NGC 5638 


16.5 


37.5 


17.0% 


NGC 7619 


40.5 


106.5 


29.7% 



- 18 - 



max=0.280345E-02 




Fig. 7.— Same as Fig. 3, but for NGC 1407. 



Fig. 8.— Same as Fig. 3, but for NGC 2986. 
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Fig. 9.— Same as Fig. 3, but for NGC 4374. 
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max=0.125187E-01 




Fig. 10.— Same as Fig. 3, but for NGC 4486. 
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Fig. 11.— Same as Fig. 3, but for NGC 5638. 
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region extends over 2 octants of the viewing sphere. 

5. Intrinsic alignment 

Intrinsic shapes of galaxies have implications for their formation and evolution, but 
the orientation of an arbitrarily chosen set of elliptical galaxies do not contain any useful 
physical information. Orientation of a very large number of ellipticals may be used to 
calculate an improved prior (Bak and Statler, 2000). In addition to above, we can also 
attempt to find correlations in orientations amongst the ellipticals of the same cluster. 

Presently, we discuss the alignment of the intrinsic principle axes of some ellipticals. It 
was assumed in Statler et al. (1999) and in Statler (2001), while investigating the intrinsic 
shapes and orientations of NGC 1700 and 3379, that the intrinsic principle axes of these 
galaxies are aligned. We now wish to examine the consistency this assumption and our 
orientation estimates. 

We choose NGC 5638, 4486 and 3379. We consider each of these in different radial bins 
and estimate the orientations of these separately. Same orientation of different bins will 
imply the alignment (Appendix A, Statler et al. 1999). The radial bins and the respective 
observational data, used in orientation estimates, are shown in table 4. 

Note that the data in bin 1 of each of these galaxies are the same as used in section 4. 
Orientations of the bin 1 of NGC 4486, 5638 and 3379 are presented in figures 10, 11 and 
13. The radial bin 2 of NGC 5638 engulfs the radial bin 1 and extends further at inner side. 
Orientation of bin 2 of NGC 5638 is presented in figure 17. We find that the plots of V of 
the orientation of bin 1 in figure 11 and of bin 2 in figure 17 are quite similar to each other. 
The peak of the orientation of bin 1 is at (#' , <f> ) = (16°. 5, 37°. 5) which is close to the peak 
(16°. 5, 40°. 5) of bin 2. These results support the notion that the intrinsic principle axes of 
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Fig. 12.— Same as Fig. 2, but for NGC 3379. 
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Fig. 13.— Same as Fig. 3, but for NGC 3379. 








Fig. 14.— Same as Fig. 3, but for NGC 4261. 
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Fig. 15.— Same as Fig. 3, but for NGC 4551. 



max=0.182070E-02 




v 

Fig. 16.— Same as Fig. 3, but for NGC 1052. 
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NGC 5638 are aligned. 

To confirm the results further, we consider bin 3, which is a narrow and overlapping 
region between bin 1 and 2. The peak of the orientation of bin 3 is (16°. 5, 37°. 5) which is 
the same as the peak orientation of bin 1. Plot of V{9\ ) of bin 3 is shown in figure 18, 
which is almost similar to plots of V of bin 1 and 2. 

A closer scruting of the plots of V reveal that while plots of bin 1 and bin 2 are very 
close to each other, plot of bin 3 is slightly different: la area of the probability of bin 3 is 
62% higher than that of bin 1. We note that Ri n is very small and R ou t is very large for bin 
1 and bin 2. However, R out of bin 3 is not very large. The values of {6 p , (/>' ) alongwith the 
area of 68% highest probability of orientations of bins 1, 2, 3 are shown in table 5. Thus, 
the photometric data and the methodology are consistent with the idea that the intrinsic 
principle axes of NGC 5638 are aligned. 

Probability V of the orientation of bin 1 and bin 2 of NGC 4486 are presented in figures 
10 and 19. The location of the maximum probability of bin 1 and bin 2 are almost the 
same. The results are consistent with the assumption of alignment of the intrinsic principle 
axes of NGC 4486. Note that both galaxies NGC 5638 and 4486, exhibit high values of |0^| 
and their orientations are well constrained. 
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Fig. 17.— Same as fig. 3, but for bin 2 of NGC 5638. 
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Table 4: Observati onal data in various radial bins used in models. 

Galaxy Bin R in R out 

tin (-out 

NGC 5638 1 8.0 36.9 0.056 0.118 8°.4 

2 4.7 36.9 0.048 0.118 10° A 

3 4.7 10.7 0.048 0.064 4°.2 
NGC 4486 1 15.0 69.0 0.033 0.100 -13°.3 

2 15.0 122.2 0.033 0.131 -15°.0 

NGC 3379 1 15.7 49.3 0.078 0.133 0°.0 

2 6.1 49.3 0.075 0.133 -4°.5 

3 6.1 15.7 0.075 0.078 -4°.7 
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Fig. 18.— Same as fig. 3, but for bin 3 of NGC 5638. 
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Table 5: Summary of the orientations of various bins of NGC 5638, 4486 and 3379. 

Galaxy Bin & cf)' p Percentage of the area of 

68% highest probability 

NGC 5638 1 16.5 37.5 17.0% 

2 16.5 40.5 13.9% 

3 16.5 37.5 27.5% 
NGC 4486 1 73.5 10.5 11.0% 

2 73.5 10.5 8.1% 

NGC 3379 1 16.5 88.5 55.4% 

2 19.5 109.5 32.3% 

3 31.5 154.5 23.7% 
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Fig. 20. — Plot of the probability distribution as a function of (#', Oy) for NGC 3379. 

Results of NGC 3379 are complex. The results of the locations of maximum probability 
and of the area of 68% highest probability of bins (1-3) are shown in table 5. We find that 
locations of maximum probability are not very close. Further, the orientation estimate of 
bin 1 is not very much likelihood dominated: la and 2a regions occupy 55.4% and 92.1%, 
respectively. 

Orientation of NGC 3379 is not well constrained even when the kinematical data is 
included (Statler, 2001, S94b). The probability distribution of NGC 3379 as a function of 9' 
and the position angle Oy of the projected y-axis of the body frame was presented in Statler 
(2001). Kinematical and photometric data were used and it was shown that while 9' is not 
well constrained, Oy is constrained. We re-examine this by using photometric models only, 
and using the data as quoted in Statler (2001). The plot is presented in figure 20, which 
can be compared with figure 5 (top) of Statler (2001). We find that the use of photometric 
data alone, worsens the situation^. Here, both 9' and Oy are not well constrained. 



1 However, the photometry alone gives satisfactory shape estimates of NGC 3379 (C08) 
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6. Results and discussion 

We have demonstrated that the photometry can be used to obtain constraints on the 
orientations of the light distribution of elliptical galaxies, especially for the galaxies with 
large position angle difference. We have obtained the possible orientations of a few elliptical 
galaxies. We have also shown that the estimates are reliable by applying our methodology 
to the synthetic galaxies, wherein the orientations are known. 

We have found that the posterior probability density of the orientation estimate is 
symmetric in pairs of octants (see, figure 2 and 12, and the relevant text). This should be 
due to the mass models. The models are reflection symmetric about the principle planes. 

We have shown that the intrinsic principle axes of NGC 5638 and 4486 are aligned by 
applying the orientation estimates using photometry. However, conclusions regarding the 
alignment of axes of NGC 3379 are not rigorously established. This is basically due to the 
reason that the orientation estimates of NGC 3379 are not likelihood dominated. 

The marginal posterior density of orientation obtained by integrating over the shape 
parameters gives the possible orientation which reproduces the photometric data most 
closely for the greatest variety of shapes. Likewise, the most probable shape obtained 
by marginalizing over the orientation reproduces the data over a greatest variety of 
orientations. These do not give simultaneous constraints of shape and orientation (Statler 
2009, private communication). 
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